De data is online beschikbaar. O.a. (hier)[https://data.overheid.nl/data/dataset/grenzen-van-alle-nederlandse-gemeenten-en-provincies/resource/d5570855-ceb7-458b-ba93-61e5ec332649]. Plaats alle shapefiles in dezelfde folder als je R script en lees de data in.
library(rgdal)
gemeentenShp <- readOGR(".","gemeenten")
## OGR data source with driver: ESRI Shapefile
## Source: ".", layer: "gemeenten"
## with 393 features
## It has 6 fields
Converteer de coordinaten naar een voor R bruikbaar format.
gemeentenShp <- spTransform(gemeentenShp, CRS("+proj=longlat +datum=WGS84"))
Zet de lijst om naar een data frame.
library(broom)
gemeentenDF <- tidy(gemeentenShp)
Bouw de kaart.
library(ggplot2)
library(ggmap)
nl<-get_map("Netherlands",zoom=8)
ggmap(nl) +
geom_polygon(aes(x=long, y=lat, group=group),
fill='grey',
size=.2,color='green',
data=gemeentenDF,
alpha=0)
Source:
D. Kahle and H. Wickham. ggmap: Spatial Visualization with ggplot2. The R Journal, 5(1), 144-161. URL http://journal.r-project.org/archive/2013-1/kahle-wickham.pdf